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FOOL. ANA PUBLICAÇÃO MENSAL, 


C. D. U. 591.756: 621.317.7 


A SHORT RESPONSE-TIME LIQUID DENSITY 
TRANSDUCER (| 


RESUMO 


Este trabalho descreve um novo anstrumento para a 
tradução continua da densidade de liquidos numa grandeza 
eléctrica (corrente continua ), cuja intensidade está relucio- 
nada linearmente com essa densidade, 

Embora seja do tipo clássico de Mutuador, este densime- 


tro é novo nas soluções encontradas para a detecção, trang- 


ferência e tradução das variações da imputado sobre o flu- 


tundor, devidas a variações da densidade do liquido em que 
aquele se encontra mergulhado. 

Daqui resultam as seguintes coracterísticas de funcio- 
namento, muito importantes num densimetro industrial : 

Pd independência das variações do coutal e da viscosi- 
dade do liquido; it) independência da pressão de trabalho; 
te) independência das variações da temperatura do liquido 
(quando desejada); ic) tempo de resposta muito pequeno; 
CA grande sensibilidade e grande reprodutibilidade; vi) 
independência da distância a que é transmitido o sinal de 
suida, 

A exactidão actual conseguida é de 0,5 *g, sendo fireil- 
mente detectáveis variações de densidade tio pequenas como 
0,0001 quem, 

Além das aplicações imediatas em todos os ramos indus - 
triais onde à medida da densidade de liquidos é necessuria 
(indústrias químicas, petroliferas, alimentares, de papel, 
de cimento, etc.) este instrumento é especiulmente mil, dao 
o seu pequeno tempo de resposta, naqueles casos em que a 
medida de densidade deve ser rápidamente obtida com um 
diminuto caudal ee liquido, É o caso, por exemplo, da me- 
didu da densidade de wm liquido obtido por dificil filtração 
de lamas de fabrico. 


bs VIRGÍLIO IGREJA 
LUIS MALTEZ 
Companhia União Fabril 


Research Center 
SYNOPSIS 
Based on a new concept of the Joater type densimeter, 


this paper introduces a new device for transducing contin- 
uousty the density of a liquid in such a way that the 


following interesting features are achicved ; 


t) independence of the Aow-rate and the viscosity of the 
liquid ; 


ti) independence 0f the working pressure ; 

vii) independence of the liquid temperature ; 

iu) a short response time ; 

v) high sensitivity and repeatability ; 

vê) enaineas in remote tranamission, 

Th: present results show an overall accuracy better than 


0,5 0 Jos. 4. 
easily detected. 


+ à deviation of MOI Sp. Er. untts bein 


Deside the úmnmediate applications in all industrial 


fields aqhere a measure of density is necessary (chemical 


and petrochemical industries, food industry, paper and 
pulp industry, ete.), the device ia very suitable in all cases 
where a continuous and quick density measurement on re- 
duced variable ow-rate sumple fa a must. Pis canhe, for 


instance, d filtrate obtained from a slurry. 


1“ Presented at the International Measurement Congress IMEKRO IV held in Warsaw, Poland, 3-8 July 1967. 
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1. INTRODUCTION 


In the present state of the art of continuously 
transducing the density of a liquid into an elec- 
trical signal, there are some features which are 
almost always desirable for the transducer. They 
are: i) independence of the liquid flow-rate and 
viscosity ; ii) independence of the working pres- 
sure; iii) independence of the liquid tempera- 
ture; iv) a short response time; v) a high sensi- 
tivity and repeatability; vi) easiness in remote 
transmission. 

The features i) to iii) are necessary to elimi- 
nate the main disturbing variables affecting the 
measure. Feature iv) is required whenever work- 
ing with very low flow-rates. A high sensiti, 
vity is a must in a number of applications and re. 
peatability is indispensable for precise measure- 
ments. At last an electrical signal capable of 
being easily remote transmitted is almost always 
desirable since the indication is generally located 
at some distance from the measuring point. 

The particular problem of continuously and 
quickly transducing into an electrical signal the 
specific gravity of a filtrate obtained from a 
slurry, at a very reduced and variable flow-rate, 
with temperature compensation and high sensi- 
tivity and repeatability, led the Research Center 
of Companhia União Fabril to develop a densi- 
meter (!) with all the above said features. This 
densimeter is based on a new concept of the 
floater type densimeter. 


(1) Portuguese patent nº 43662 of the 2.26 1965, in the 
name of Companhia União Fabril. 
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In order to obtain independence of the liquid 
flow-rate and viscosity, the floater is located 
along the current lines, these being horizontal; 
in this way the disturbing forces do not affect 
the floater apparent weight, vertically measured, 
which is related with the liquid density. The 
temperature independence is obtained through a 
reliable hydraulic system, installed inside the 
floater. The use of a tubular floater, concentric 
with a closed tube conducting the fluid, both 
parts being of almost equal diameters, allows the 
measuring in closed pipe on a reduced liquid 
volume ; consequently a short response time and 
independence of the liquid pressure are obtained. 
A force balance weighing system wich does not 
mechanically touch the floater (the force trans- 
mission is carried out magnetically) guarantees 
a high repeatability. Finally, the output signal 
(0-10 mA, D.C.) of this force balance weighing 
system, being the restoring variable through its 
feedback circuit, may be easily remote transmit- 
ted without affecting the measure. 


2. WORKING PRINCIPLE 
2.1. Densimeter 


Asshown in figure 1 this densimeter is basic- 
ally composed of two units mechanically con- 
nected. The first one is composed of a tubular 
floater f concentric with a nonferromagnetic 
material tube € coupled through an articulation 
O. The tube C is the measuring vessel and it is 
meant to conduct the liquid. The floater thor- 
oughly immerges into the fluid. 
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The other unit is composed of an automatically 
balanced weighing set having a lever t, to which 
a ferromagnetic core ni, a permanent-magnet n> 
and also a permanent-magnet or a permanent- 
-magnet set ms are attached. Besides, the weighing 
assembly also disposes of two coils Bi; and 
Bs, the forme: being the detection coil of the ni 
core position and the latter the driving coil which 
exerts an electromagnetic force on the permanent- 
-magnet ns. 

These coils are rigidly linked to the tube C by 
the parts pointed out in figure 1 as |, and 1,. 

Inside the floater £ lays a permanent-magnet 
m, in such a position that it applies a repulsive 
force on the above said permanent-magnet set 
m,. The screw p is the limiter of the movement 
of the weighing set lever and f, is its fulcrum, 
The coils are electrically connected according to 
the diagram of figure 2 where A represents a 
regulating device (an electronic amplifier in its 
simpler form); mA represents an ammeter. 


Fr 
1 
] 
| 
l 


a. 


magnetic repulsion force R due to the magnets 
m, and m,, upwards. 
From the torque balancing we have: 


P.oa=l.b+-R.x (1) 
where a, b and c are the arms corresponding to 
each force and 

[== Vu (2) 
V being the floater volume and u the liquid spe- 
cific gravity (weight of the volume unit). When- 
ever the liquid specific gravity changes the floater 


slightly moves and the repulsive force R chan- 
ges too. From (1) and (2) we get 


AR =— K,.du (3) 
where AR is the change of the repulsive force 


R, du the corresponding change of the specific 
gravity u and K, a constant obtained from the 


t 


Fig. 2 


The working principle of this densimeter is, 
in short, the continuous weighing of the floater 
immersed in the liquid by means of an automat- 
ically balanced electromagnetic weighing set 
which does not mechanically contact the floater. 

As the apparent weight, which unbalances the 
weighing set, depends on the liquid specific grav- 
ity, the intensity of the rebalancing current may 
be related to the specific gravity of the liquid 
where the floater immerges. This current may be 
used as a measure for the specific gravity of the 
liquid. 

As a matter of fact, the floater is balanced 
under the action of three forces: self weight P 
dowwnards, impulsion I due to the liquid, and 
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constants of (1) and (2). As R changes the 
weighing set endures a small unbalance which 
causes a small displacement of the coren, inside 
the coil B,. In this coil, which is a differential 
transformer, the displacement gives place to an 
electrical voltage e with the frequency of the 
voltage E which feeds the primary winding of 
the mentioned differential transformer. Voltage 
e, by its turn, is the imput of an electronic am- 
plifer A which sends a continuous corrent i to 
the driving coil B,. This current gives rise to a 
force f, on the permanent-magnet m, and, there- 
fore, on the weighing set lever. This lever suf- 
fers the action of three forces: self weight P”, 
downwards (at the fulcrum left or right according 
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to the position of its center of gravity), repulsive 
force R and driving coil electromagnetic force f,. 

From the torque balancing equation and fol- 
lowing a reasoning similar to the one that led us 
to equation (3), we can draw the following con- 
clusion : 


BE == —K, dB (4) 


where Ks is a constant. 
From (3) and (4] we have 


AL =K,.K,.du (5) 


i. e., the variations of the force f, are propor- 
tional to the specific gravity variations of the 
liquid where the floater is sunk. On the other 
hand as 


Ai=K,.4f, (6) 


where K, is a constant, is the equation that re- 
lates the d.c. intensity i, feeding the coil B,, with 
the force f, applied to the permanent-magnet n,, 
we obtain from (5) and (6) the equation 


Ai=K .K. Kd (7) 


which is the basic equation of the densimeter. 
2.2. Temperature compensation system 


A well known method for temperature com- 
pensation is the use of a temperature transducer, 
the output of which is added to the densimeter 
one at the work temperature, in such a way that 
the final signal corresponds to the liquid specific 
gravity at the reference temperature. For very 
accurate measurements, this method implies a 
compensation system of high quality which 
increases the price and complexity of the device. 
This complexity generally implies a more expen- 
sive maintenance and a certain loss of reliability. 
On the other hand the feedback current of the 
weighing set related with the liquid specific grav- 
ity should preferably be completely free from 
any disturbing signals. For this purpose a par- 
ticular temperatuce compensation system (*) 
well suitable to this densimeter was developed 


(*) Portuguese patent nº 46 649 of the 10,29,1966, in 
the name of Companhia União Fabril. 
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wich renders the device insensible to those fluc- 
tuations of specific gravity exclusively deriving 
from changes in the liquid temperature. 

The working principle of this system is to 
obtain a correspondence between the position 
of a mass placed inside the floater and the tem- 
perature of the liquid, in such a way that the 
changes in the buoyancy torque, due to fluctua- 
tions of the liquid specific gravity (caused by the 
changes in its temperature), are exactly compen- 
sated by the changes in the floater weight torque 
obtained from the adequate displacement of that 
aforesaid mass. In fact, from the above equations 


Poa=1I.b+R.c (1) 
Il=V.u (2) 

we get 
P.a=Vu.b+R.c (8) 


where P is the floater actual weigth, Vis the 
floater volume, u is the specific gravity of the 
liquid, R is the magnetic repulsion acting on 
the perinanent magnet fixed to the floater and a, 
b and c are the arms corresponding to each force 
and referred to the axis O. In order that R does 
not change with the fluctuations of the specific 
gravity (caused by changes in the liquid temper- 
ature) and, in this manner, does not change the 
feedback current of the self-balanced weighing 
set (that is to say, does not change the output 
of the device), the floater center of gravity must 
suffer the displacement da given by the equa- 
tion : 

P.da= V.b. dy (9) 


On the other hand, we have 
je SO dA T (10) 


where K is the variation coefficient of density 
with temperature. 
From (9) and (10) we get: 


Pebg=V.b.&K,aT (11) 


In brief, we can say that the device output 
signal does not change with the temperature 
fluctuations AT if, simultaneously, the center of 
gravity of the floater is displaced of da accord- 
ing to the equation (11), which is the basic 
equation of the temperature compensation 
arrangement. 
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Fig. 3 


The mass placed inside the floater is a liquid 
which fills up a system formed by two chambers 
linked by a tube of very small diameter (fig. 3). 
One of the chambers has a constant volume 
(rigid walls); the other one has a variable vol- 
ume (movable walls). The device is fixed inside 
the floater, the variable volume chamber being 
at the side of the floater pivot and the constant 
volume chamber at the opposite side. When the 
surrounding temperature changes the quantity 
AT, the interior liquid volume must change the 
quantity Av, in such a way that 


PA a l,dv u (12) 
where | is the distance between the chambers 
(approximately constant) and u” is the final spe- 
cific gravity of the liquid filling up the system. 
As we have, on the other hand: 


Av=kK'.4T.v, (13) 


where v, is the initial volume of the interior 
liquid and K” its relative variation coefficient of 
me with temperature, we get from (11), (12) 
and (13): 


Lu.K .y=V.b.K (14) 


This equation shows that it is possible to 
obtain a total temperature compensation ad- 
justing v, or |, when uí and K' are fixed by the 
choice of the interior liquid of the system. 

Fig. 4 is a block diagram illustrating the 
feedback circuits in the device. 


 P | 
|) 
uu a4 x | 
A . 
T | floater [O 
sm | ) 
| | 


xo : 
R ' 
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3. DETAILS OF OPERATION 
3.1, Range and Zero Adjustment 


The instrument output is a D. €C. current 
varving from zero to 10 mA. The sensitivity is 
primarily fixed by the floater volume V; adjust- 
ments are obtained through a «range resistance» 
shunting the driving coil. Zero adjustments can 
be made changing the position of the weighing 
set counterweight CC... 


3.2. Independence of Flow Rate and Viscosity 


A complete independence of flow rate and 
viscosity is obtained only when the floater is 
perfectly concentric with the pipe conducting 
the liquid, However if the repulsion force R 
between magnets m, and m, is large enough, 
the floater position is practically fixed for all 
values of density in the range and so a very 
satisfactory independence is attained. 


3.3, Calibration of the Temperature Com- 
pensation System 


As shown at the end of section 2.2 this 
adjustment can be made by changing the values 
of v, (the volume of the compensation liquid) or 
| (the distance between the chambers). Both 
methods have already been used the latter one 
permitting an easier calibration. In this method 
a metallic flexible tube has been used. 
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DENSIMETER WORKING PRINCIPLE 


! Jo (do) 


(d-VARIABLE , T-CONSTANT) 


1) d-CONSTANT 
d=d, —— I=Jo-—— R=Po — F=b — Í=lo 


bati cr: | 


EE ad E md Sa 
á » 
" % 
E a ba) 
a 
É x 3 
“e ei nd E 
4 é 
Dc e Sp e IC aa st e ai 
r + 
E ] ê E 
maca Nm a = E e BR | + 
| | 
Leal 


d=do=CONST = 


É x Eat EA 
L 
4 


d=do = COMST 


mas == — — — — e e mm 


qr) = 


> a() ————————— 


Ê | dois | 


— 2) d-INCREASES | 
[a a d>d, —— [ >» ——P<R, ——/ >fs Eá (>lo 
d>do emma i>lo 


e d ——s» | 
[feio ea 
Ie — 3)d-DECREASES 


dedo —— I<lo—— R>R, —— [<f— i<lo 
d< do a o t<i, 


DENSIMETER WORKING PRINCIPLE 
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In high sensitivity measurements the use of a 
temperature compensation system is inevitable 
due to the slight influence of temperature on 
the floater dimensions. 


3.4. Limitations 


Any disturbance affecting the floater apparent 
weight, alters the output. Expected disturbances 
are generally gaseous bubbles and solid particles 
which adhere to the floater. Nevertheless the pipe 
and floater geometry permits to work at a liquid 
velocity which reduces that adherence. Anyway 
it is advisable to avoid such a problem. 

Another important disturbance may be the 
floater corrosion because this affects the floater 
weight and/or volume. This generally leads to a 
zero shift. The floater material must therefore be 
carefully chosen according to the liquid. 

Friction torques at the floater articulation are 
negligible compared with the liquid impulsion 
torque variation. 


IO | 


output current 
Ea 


0. 7200 
specific gravity 


4. TEST RESULTS 


Diagram I shows the results obtained from 
different mixtures of aircraft fuels using an in- 
strument without temperature compensation. The 
actual temperature was T=21+1ºC. It was 
possible to apreciate specific gravity changes of 
0.0003 g cm” *in the range umin = 0.7059 g cm” *, 
Umax == 0.7418 g cm”? (du = 0,.0359 g cm” *). The 
twenty measured values are located between 15 
and 80'/, f.s.d., approximately. The resulting 
maximum error was 1.35º/, f.s.d. referred to the 
best straigth line adjusted to the points. The 
floater used had a lenght of 350 mm and a di- 
ameter of 26 mm, the inner diameter of the pipe 
conducting the liquid being 34 mm. With a pul- 
sating flow-rate of 4 | min”! furnished by a 
peristaltic pump, which did not disturb the out- 
put, those dimensions led to a response-time of 
about 5 seconds. 


0. 7300 0. 7400 
u tg. cm?) 


Diagrama 1 
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Hig. 7 — Primary element 


Fig. 8 — Indicator, electronic set and primary element 


5. CONCLUSIONS 


The interesting features of this instrument are 
its simplicity and the quality of its output signal 
even in presence of strong disturbances, as for 
instance, pulsating flow. This results mainly from 
the particular geometry and position of the floater 
and from the fact that the density is transduced, 
without mechanical contact, into a force a varia- 
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ble which can then be easily transduced into an 
electrical or pneumatic signal, using well-known 
feedback techniques. The soundness of its work- 
king principle has been confirmed through the 
satisfactory behaviour of several instruments 
already in service, even of those which were not 
very carefully construted. It is therefore reason- 
able to foresee a promising future for this de- 
vice. 
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DA TEORIA DA FONTE MÓVEL NUM MEIO INFINITO 
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RESUMO SYNOPSIS 

Utilizando a transformada de Laplace, generaliza-se, A solution for transient regime temperature distribu- 
para o regime transiente, q expressão da distribuição de tion in an infinite solid heated by a constant velocity 
temperaturas originada num corpo infinito por uma fonte moving heat source is found by means of Laprace trans- 
de calor móvel com velocidade constante, form, 


1 — Introdução 


Resumem-se, em linhas gerais, aspectos do regime estacionário do fenómeno e apresenta-se, 
como exemplo de aplicação da transformada de Laplace — ferramenta das mais poderosas na análise 
teórica dos aspectos transientes dos fenómenos de condução de calor —, uma extensão ao estudo 
do regime transiente. 


2— Modelo teórico 


Considera-se uma fonte pontual de calor (libertando uma quantidade de calor Q, na unidade 
de tempo, para t >> 0) que se desloca a uma velocidade V constante (pressuposta de módulo bas- 
tante inferior ao da velocidade da luz), para t > o, um meio isótropo, infinito em todos os sentidos. 
Admite-se, mais, que o calor introduzido no meio, lhe não altera as propriedades térmicas, não ori- 
ginando, também, qualquer mudança de fase. 

Associa-se ao meio um referencial fixo (O,x,y,z) é à fonte um referencial móvel (Op xr yR ZR) 
tais que 


para E=0 Or=O Ox = Oxp Oy = Oyr Uz = LVDzr 
para t>o Ox = Oxp Oy // Oyr Oz //Ozr 


2.1 — Tratamento matemático 


O tratamento clássico, que, aliás, adoptaremos, visa a determinação do campo de temperaturas 
associado ao referencial móvel, o que exige a integração da equação do calor escrita em coordenadas 
deste referencial. 

Por demais conhecida é a técnica de passagem do referencial fixo ao móvel [1], para veloci- 
dades muito inferiores à da luz, pelo que nos abstemos de a transcrever, partindo já da equação 


conveniente. 
dº T d* T d* T V dodT 1 0 E 


jo ao dp qa “fe = (1) 
dd XR J YR TR) Zp zu XR x U t 
x — difusidade térmica t — tempo T— temperatura 


— Passaremos a substituir, xr, YR, ZR Por x, y, Z, respectivamente, uma vez que, beneficiando 
de uma maior simplicidade de escrita, se não prejudica o texto com qualquer ambiguidade. 


2.1.1 — Regime estacionário 


Postula-se, na análise deste, que suficiente tempo decorreu já a partir do instante de arranque 
(t==0t), para que o campo de temperaturas tenha atingido o seu perfil de estacionaridade. Tal 


suposição equivale a admitir que se atingiu uma altura em que, em cada ponto (x, y, z), as flutuações 
— E n * É jà 
de T com t são desprezáveis, ou seja, que — = 0, 
ut 
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A equação cuja integração se exige sera, então 


9º T 2 T RE Va 
ima (3) 
d x? d y* Jd z a | 


Tratou, originalmente, este caso RosextHaL [1). A solução obtida por este autor, através de 
um processo de integração artificioso, mas legitimo, reveste a forma 


V 
-— (x +) 
27 
E (3) 
arkr 
k — condutividade térmica T-—temp.º (origem: o graus) 


p= xyz 


A contracção do conjunto de variáveis x, y, z para x, r justifica-se pela simetria axial do 
fenómeno. 


2.1.2 — Regime transiente 


2.1.2.1. — Generalidades 


No que se disse até agora, uma lacuna subsiste: a maneira como evolui o campo de temperatu- 
ras, desde o instante inicial (temperatura nula uniforme) até ao regime dito estacionário atrás descrito. 


7 ) 
Pretende-se conhecer processos de «deformação» das isotérmicas no domínio de tempo em que ad 
apresenta, em cada ponto, valores não desprezáveis. Haverá, então, que integrar 
PI AT. L,VoT IT 
: + | ns >= —— (4) x 
od xº dy? d z* x dx a dt 
2.1.2.2 — Integração de (4) 
1) Estabelecimento das condições fronteiras (no espaço e no tempo) 
= 0 T T=0 em todo o domínio 
t qualquer T+ O quando | * oo 
r-— 
REA E 
2) Proposição da primitiva T=e * firt)=e * fxy,zst) 
(4) degenera em 
)? E o* E a A RR 1 0£ 
Cp +- (0) fe (5) 
d x* o ye dz” 2% | 4 dt 
3) Aplicação da transformada de Laplace em ordem ao tempo. 
Tomando 
E oa — st 
t (x,y,2,5) = À f (x,y,z;t) e dt=[, 
E 
(5) conduz a 
o E d? E 2 f | Vy- E Dara 
e je came fl — — 60 E E iss — T(o*,x,w;z) (6) 
d x* dy* d zº 2% x a 


E =" + * . = | n o 
(*) Zinsmeister e Dixon estudaram, pormenorizadamente, o problema, no caso do corpo infinito a uma dimen- 
são, aplicando, também, a transformação de Laplace /2). 
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4) Introdução da condição inicial 
T(ot,x,y,2) =0 


92 E o? E da É Vi? s 1- 
—— + >|(—) += |f=0 
o x* E dy? dz (5) | (0) 


5) Mudança de coordenadas E (x,y,2,8) — f(r,s) 


(7) transforma-se em 


df 2. df | ( V y s ]- 
es são Re: ques =— + — | f=0 8 
dr r dr | 2%) E á 
6) Proposição da primitiva E(r,s) == 8 (1,8) 
r 
(8) origina 
d”, VA | 
(SJ + Leme o 
dr” [122 - 
7) Integração de (9) 
Ss (VN? (Via 
É — dm r cá so PR: 
VE) VS a (10) 
g = Cie + Cse 
8) Introdução da segunda condição fronteira. 
v 
ao * 
Posto que r é um infinitamente grande de ordem inferior a e si e deverá ser 


lim T=0, somos levados a concluir que C/==0 para que a condição seja satisfeita. Teremos, por- 


p—+ 00 


tanto 


O 


IS a (NY 
dia Ei a (11) 


9) Percorrendo, em sentido inverso, as alíneas anteriores e desfazendo as transformações da 
forma da função que nelas se operaram, chegamos a 


aa ga = o 
a A a o (12) 


10) Quando V=0, o regime há-de degenerar na solução para uma fonte parada, libertando 
a quantidade de calor Q por unidade de tempo 


É Q GEE = == (13) 
axkr V 4at 
que tem por transformada 
+ 
-— | Eu I 

T= E e Q 

s 4mkr 

Então, será nodes 
=/ Es -V E 
1 z Q 1 a 
—€ = .. 

r akrr 5 
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11) cri Ar hdi (14) 


12) Determinação da transformada inversa, em que está, apenas, interessado o grupo de 


factores 


| V q 
ex —— p lerfe | -—= +— ani |2 
Li | 2V at 2 a | [2] 
13) 
V V V 
e É 1 is ES V AR Za" r V VE 
T= q e le erte| Fa sã f*|+te erfc|) —== + =— A 
akz r 2Vat 2 ) x 2V at 2V x 
(15) 


2.2 — Notas finais 


1) A integração de (4) tem fortes pontos de contacto com a de (2), mostrando-se ainda, vali- 
dos muitos dos artifícios usados nesta. 

2) (15) confirma as condições fronteiras impostas e satisfaz à equação às derivadas parciais 
de partida. Além disso, degenera na expressão do regime estacionário quando t — cs. 

3) (3) e (15) encontram aplicação imediata em problemas de carácter tecnológico, dos quais se 
salienta a soldadura por arco eléctrico, fornecendo um modelo teórico razoavelmente aceitável para 
o estudo da distribuição de temperaturas. À isto nos referiremos noutro artigo. 
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REPRESENTAÇÃO TENSORIAL DO COMPORTAMENTO 
INCREMENTAL DE ELEMENTOS DE REDES ELÉCTRICAS. 
ALGUMAS APLICAÇÕES (*) 


RESUMO 


Neste artigo introduz-se a ideia de definir por um mé- 
todo tensorial o comportamento incremental de elementos de 
redes eléctricas, 

Sião definidos os tensores impedância e admitância im- 
cremental e são estabelecidas as relações fundamentais num 
circuito de elementos cujo comportamento é definido por ten- 
sores, 

Indicam-se algumas aplicações dos métodos expostos na 
definição do comportamento incremental dos circuitos e apre- 
senta-se um exemplo de aplicação à análise de variações 
de tensões, 

Os conceitos de tensor impedância e adimitância incre- 
mental mostraram-se muito cômodos e úteis, por exemplo em 
problemas de variações de tensão, variações de cargas, per- 
das e de comportamento de redes tendo em conta a actuação 
dos elementos de regulação, 


o. INTRODUÇÃO 


Neste artigo introduz-se a ideia de definir por 
um processo tensorial o comportamento incre- 
mental de elementos de redes eléctricas. 

Há elevado número de problemas concretos no 
estudo de redes, designadamente tendo em conta 
a actuação dos sistemas de regulação, em que é 
necessário ou útil definir o comportamento da 
rede na vizinhança de uma situação conhecida, 
para analisar o efeito de uma pequena alteração 
de situação. 

É em muitos casos inaplicável a hipótese de 
linearidade dos elementos do sistema. 

Verificou-se ser em muitos casos possível e 
util fazer uma análise do comportamento incre- 


por CARLOS DE MEDEIROS PORTELA 
Doutor em Engenharia Electrotécnica 
Chefe do Centro de Cálculo da Companhia 
Nacional de Electricidade 


SYNOPSIS 


In this paper is introduced the idea of defining by a ten- 
sorial method the incremental behaviour ofelectric networks 
elements, 

The impedance and admittance incremental tensors are 
defined, and are established the fundamental! relations in a 
circuit of elements the behaviour of which is defined by ten- 
sora. 

Some applications are indicated when defining the cir- 
cuits incremental behaviour and an example of the appli- 
cation to the voltage variation analysis is shown. 

The concepts of impedance and admittance incremental 
tensors were found very convenient and useful, for example 
on the problems of voltage variations, loads variations, los- 
ses, and networks behaviour considering the actuation of 
regulation elements, 


mental considerando um primeiro termo linear 
desse comportamento, a sobrepor ao estado base 
anterior, hipótese bastante mais correcta que a 
suposição de linearidade do sistema ou outras 
hipóteses simplificativas correntes. 

Na análise deste comportamento incremental 
não é em geral aplicável o conceito de impedân- 
cia (ou admitância), como cociente entre os com- 
plexos representativos do acréscimo de tensão e 
do acréscimo de corrente, por essa impedância 
ser intrinsecamente função da tensão incremental 
(ou da corrente, ou simultâneamente da tensão 
e corrente). No entanto, é em geral aplicável o 
conceito de tensor impedância (ou admitância) 
incremental, definidor da relação entre os vec- 
tores tensão e corrente incrementais. 


(*) A matéria deste artigo foi objecto de uma comunicação ao Symposion sur I' Analyse et la Synthêse des 
Réseaux Electriques, realizado em Bucareste em Outubro de 1967, no qual o autor tomou parte a convite da Academia 


das Ciências da Roménia. 
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1. OS TENSORES IMPEDÂNCIA E ADMI- 
TÂNCIA INCREMENTAIS 


1.1 — Definições (!) 


Seja um dipolo, em regime alternado sinusoi- 
dal, caracterizado numa determinada situação de 
referência por uma tensão U, e uma corrente Ir, 
sinusoidais, representadas, num determinado re- 
ferencial, pelos correspondentes vectores, U e |, 
a duas dimensões a e b, de componentes Us, Up 
e la, lb. 

Consideremos uma situação vizinha deste di- 
polo caracterizada pela tensão U=U, + AU e 
pela corrente [=1, + Al, sendo SU e 4] «pe- 
quenos». 

Defina-se o tensor impedância incremental, Z, 
por 

ôU = Z,41 


ou 
AUa = Zaa . Ala + Zab . Ap 


AUp = Zpa . Ala + Zppb . Alb, 


a menos termos de ordem superior. 
Defina-se, análogamente, o tensor admitância 
incremental, Y, por 


8] =Y. 4U, 


ou 


Ala ma Yaa . àUa + Yab e AUb 


Alb = Yba . AU, je Ybb . Us, 
a menos termos de ordem superior. 

Este tipo de representação engloba como caso 
particular a do comportamento incremental das 
máquinas rotativas, considerando comportamentos 
diferentes para efeito de grandezas longitudinais 
e transversais, e permite ter facilmente em conta 
diversas formas de regulação. 

A título de exemplo, indica-se a seguir o valor 
destes tensores de comportamento incremental 
para dois tipos de comportamento de circuitos. 


(!) Por simplicidade de notação, não se explicita o 
carácter covariante ou contravariante dos tensores utili- 
zados. É imediato, no entanto, o estabelecimento de um 
conjunto coerente de definições de covariância e contra- 
variância, 


84 


1.2 — Tensores impedância incremental e 
admitância incremental em circuitos 
em que são definidas elasticidades da 
potência activa e da potência reactiva 
em relação à tensão 


Em muitos elementos de circuitos, são conhe- 
cidas, exacta ou aproximadamente, as elastici- 
dades da potência activa e da potência reac- 
tiva em relação à tensão, i. e., na vizinhança de 
um estado tem-se (!) 


sendo Pe Q a potência activa e reactiva absor- 
vidas pelo dipolo a cujos terminais está aplicada 
a tensão: 


P = Uala + Ublp 
Q = Ubpla — Ualb . 


Sob outra forma, numa vizinhança suficiente- 
mente pequena da situação base, de referência, 
referenciada pelo índice r, a menos termos de 
ordem superior, tem-se 


P JU| Ns 
—— = — — 1+-: 
Cu) 


[UI — (Ulr 
Ur á 


Designadamente, no caso de cargas correspon- 
dentes ao conjunto de uma rede de distribuicão, 
é em geral fácil obter, por determinações expe- 
ditas (variando a tensão aplicada), os valores de 
: en. Conhecendo o tipo de rede (proporção 
aproximada de carga de iluminação, industrial, 
etc.), é também em geral possível uma estimativa 
dos valores de : e 1. 

Como é evidente, no caso particular de impe- 
dâncias lineares passivas, tem-se 


a elasticidade de uma 
EY 
variável genérica X em relação a uma variável genérica Y. 
Representa-se por |X| o módolo de um vector genérico X, 


(1) Representa-se por 


X| = [x + x 
V Ca b 
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Em. 
Sendo 
P=P, + à P 


Q = Or + Á t, 

Pr = Uala + Ublb 

O, = Usla — Ualh 
AP = Ua Ala + Up4l, + 14404 + 1,4Up 
3Q = Updla — Uadlp + 144Up — 1,dU.. 


Fazendo os cálculos, obtém-se, a menos termos 
de ordem superior: 
1 


Pp (E—1) + Q!(n—1)' 


Zaa 
- | Uala [uí + U. (1+:—n)] + 


+ Uplo [UÉ (11) + UC (6—2)] | 


| -U, a[UZ—1) + Uj(n— 1] + 
+U, [US+ US e+n)| 
E ua 1 


Pe -D+AHO(n—1) 


|U, [UG =: +m) + UM 


— U,1, [Us(s—1) + U? ( -1)| 


1 


Z a E == Eee 
bb p? (e—-1) + Q'(n—1) 
: | U, la UU, (E—I)+ U. (n— 1)! + 
+ Up 1, [UÉ (1 +:—n)+ US] | 
1 | | ” 
=—s E 1 FEM fr so t(n—1)] 
Em ul: [Ou a LUZ ( 1) + U, ( 1) + 
+U, 1 [US (1 + em) + UZ) 
1 | ] | 
Yap e [0 /ER “| U, Ea [UÍ (E—1) + Uf (1 1) ) — 


— U, Ip [UZ + UZ (1 =" +) 
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Ya = [= Uol [UR =) + URIA 
+U, 1, [UL(—1) +UZ 6-1] 
Yop = E (Us [US +Uga+ en] + 


+U, 1 [UZ (11) + U, 610] | 


No caso particular de um sistema de eixos p, q, 
tal que p coincida com a direcção da tensão em 
condições base, vem: 


Y tp (e — 3) 
Pp U sAsiré me 


p 
Hg E = E 
ap (ma) 
» 
Cag E 


Isto é: 


a) Para efeitos de tensão incremental perpen- 
dicular à tensão base, o tensor tem o mes- 
mo efeito da admitância definida pelo co- 
ciente da corrente e tensão em condições 
base, em representação complexa. 

b) Para efeitos de tensão incremental paralela 
à tensão base: 


b.1) a admitância incremental na direcção 
da tensão base é igual ao valor corres- 
pondente à relação entre corrente e 
tensão em condições base multiplicada 
por (c-1); 


b.2) a admitância incremental na direcção 
ortogonal à da tensão base é igual ao 
valor correspondente à relação entre 
corrente e tensão em condições base 
multiplicada por (1-1). 


A condição a) traduz o facto de, para efeitos 
de rotações da tensão, haver uma rotação igual 
da corrente sem variação de amplitude. 

A condição b) traduz o facto de, para variação 
de amplitude da tensão, sem mudança de direc- 
ção, a componente activa da corrente ter uma 
variação relativa igual à variação relativa da ten- 
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são multiplicada por (:-1) e a componente reac- 
tiva da corrente ter uma variação relativa igual à 
variação relativa da tensão multiplicada por (1-1). 


1.3 — Tensor impedância incremental de um 
circuito caracterizado por tensão de 
amplitude constante e potência activa 
constante 


No caso de um gerador com regulação ideal 
de tensão e potência activa, o módulo da tensão 
e a potência activa são constantes. 


Fazendo os cálculos obtêm-se: 


U, U, 
aa Q 
U, 
Lab FR SA 
Q 
U* 
dei Q 
U, U, 
Lp E mg 
Q 


O determinante deste tensor impedância é nulo 
e não é possível definir um tensor admitância. 
Este facto resulta de, num circuito nestas con- 
dições, haver uma infinidade de hipóteses de va- 
riação de corrente que correspondem à mesma 
variação de tensão; a corrente define a tensão, mas 
a tensão não define a corrente. 


hd 


RELAÇÕES FUNDAMENTAIS NUM CIR- 
CUITO DE ELEMENTOS CUJO COM- 
PORTAMENTO É DEFINIDO POR TEN- 
SORES 


Consideremos um circuito, constituído por 
vários elementos, cujo comportamento em regime 
alternado sinusoidal seja traduzido por tensores 


.. A . 2; À = . 
impedância ou admitância incrementale por f.e.m,, 


ligando os nós deste circuito, como representado 
na fig. 1, e consideremos correntes de circulação 
nas malhas, le. 

Considerando para as f.e.m. e correntes as 
diferenças entre dois estados vizinhos, sendo AE; 
a diferença entre as somas da f.e.m. na malha k 
nesses dois estados, tem-se: 
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Zu Al =++ Zu Alo + a + Zim Am ue SE, 


ZM St + Z32 Bla e o o a Zam Am = DEs 
Zml dl Zm? Als + Es A 5 +- Zm A Im = AEm, 
Es E, + E. 


Fig. 1 


sendo 


Zi o tensor impedância incremental pro- 
pria da malha i, soma dos tensores 
impedância incremental ao longo da ma- 


lha i, 


Zii o tensor impedância incremental mutua 
entre as malhas i e j, soma dos tensores 
impedância incremental nos troços co- 
muns às malhas i e j, multiplicados por — 1 
no caso de Il; e 1, terem sentidos con- 
trários nesses troços. 


Tem-se assim um sistema de m equações vec- 
toriais a duas componentes, de coeficientes ten- 
soriais, ou de 2 m equações escalares de coefi- 
cientes reais, sendo m o número de malhas, sis- 
tema de equações que, para uma variação de 
f.e.m. definida por AE, define as Variações 
de corrente nas diferentes malhas, Ali. Conhe- 
cidos os valores de 4lh, pode obter-se directa- 
mente qualquer tensão ou corrente incremental- 

Ter-se-á 

Zi =— Zi, 


sendo portanto simétrico o sistema de equações 
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de coeficientes tensoriais, embora não sendo em 
geral simétrico cada um dos coeficientes Z, |; con- 
sequentemente, o sistema de 2 m equações esca- 
lares não é simétrico. 

É possível generalizar os conceitos anteriores 
e considerar uma rede com elementos activos, em 
que o valor de uma grandeza da rede (por exem- 
plo tensão, corrente ou combinação linear de ten- 
sões e correntes) vá influenciar de maneira não 
simétrica uma ou mais tensões ou correntes ou 
combinações lineares (!). A consideração de ele- 
mentos deste tipo traduz-se, no sistema de equa- 
ções, pela não verificação, com generalidade, da 
igualdade 

Li, = Lj 
Aparte isso, no entanto, mantém-se o aspecto 
formal do sistema de equações, com alterações 
do significado dos coeficientes. 

Consideremos agora uma formulação dual das 
relações do circuito, sendo U, as tensões dos 
nós k do sistema em relação a um nó de refe- 
rência (0), sendo Y, o tensor admitância incre- 
mental das ligações entre os nós i e j, sendo Y,; 
a soma das admitâncias de todas as ligações en- 
tre o nó ie os restantes nós incluindo o nó 0, 
de referência, e sendo À Ik a variação de corrente 
exterior injectada no nó k (fig.2). 


Fig. 2 
Tem-se 
Yu AU, ii Yj AU, — Ed Eis AU, = É) , 
— Yy AU, + Ya AU, — ...— Ya, AU, =, 
= Ya AU, Es Tas AU, — eee TF à (PER AU, eee dl, 


(!) Para efeitos incrementais, a potência activa e reac- 
tiva e grandezas semelhantes podem ser consideradas 
combinações lineares de tensões e correntes ou das res- 
pectivas componentes. 
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ou seja, um sistema de n equações vectoriais de 
coeficientes tensoriais ou um sistema de 2 n equa- 
ções escalares de coeficientes reais. 

No caso de elementos passivos, será 


Yy = Yi 

A consideração de elementos activos, actuando 

de maneira não simétrica, implica que esta úl- 

tima condição não seja satisfeita, mantendo-se no 

entanto o aspecto formal do sistema de equações, 
com alteração do significado dos coeficientes. 


3. APLICAÇÃO DOS MÉTODOS TENSORIAIS 
DE DEFINIÇÃO DE COMPORTAMENTO 
INCREMENTAL 


Elaborou-se um conjunto de programas para 
ordenador digital que permitem: 


a) Definir os tensores admitância e impedância 
incremental para várias hipóteses de com- 
portamento dos circuitos, a partir das con- 
dições base. 

b) Determinar as grandezas num circuito de 
elementos definidos por tensores admitância 
(ou impedância), por um métódo iterativo. 

c) Reduzir um circuito de n nós a um circuito 
de menos nós. 

d) Determinar as grandezas por um processo 
de eliminação. 

e) Obter a matriz inversa do sistema de 
equações de admitâncias ou impedâncias 
e consequentemente determinar : 


— amatriz das admitâncias detransferência Yi 
entre a f.e.m. na malha k e a corrente na 
malha j; 

— a matriz das impedâncias de transferência . 
entre a corrente no nó k e a tensão no nó j. 

f) Determinar as grandezas que interessam 
para diversas aplicações, por exemplo, os 
acréscimos de perdas e as variações de ten- 
são (referidas às tensões base), para uma 
determinada variação de correntes ou ten- 
sões aplicadas a alguns nós. 


Estes métodos e programas podem, ôbviamente, 
ser utilizados para representar circuitos com ele- 
mentos representados por tensores (independen- 
temente terem ou não o significado incremental) 
e, em particular, para representar elementos cor- 
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respondentes a impedâncias no sentido normal, 
caso em que 


Lpp = Las 
Lra = — Lab 
Top Lua 
Yia=—Y 


A utilização destes métodos e programas mos- 
trou-se útil e cómoda no estudo do comporta- 
mento de redes, designadamente em problemas 
de perdas e de perturbações e variações de tensão. 

Verificou-se, por outro lado, que os resultados 
obtidos em alguns casos são apreciavelmente 
diferentes dos resultados obtidos com alguns 
métodos correntes, por exemplo representando 
os circuitos por impedâncias constantes, sendo 
portanto necessário, nesses casos, a adopção de 
métodos diferentes, do tipo dos descritos aqui. 

Para o cálculo iterativo, considerou-se como 
base, um processo de iterações nodais, do tipo 
Gauss Seidel. 

Considere-se a equação 


— Yj dU, — Y,; AU, —... + 
+Y, àU, — ads a SU, =d1, 
donde, sendo Y,—! o inverso de Y,., 


/ n 
il 
+ | 


Numa dada iteração, é calculado, por meio da 
relação anterior, o valor AU; a partir dos va- 
lores de ÀU nos restantes nós. Sendo AU; o valor 
anterior de AU;,, é dado a ÁU; um novo valor, 


SU, — AU, + x, (3U;, Ea AU; o p 


sendo « a um factor de aceleração de conver- 
gência (nalguns casos substituído por um tensor 
de segunda ordem). 

Na maior parte dos casos, escolhendo experi- 
mentalmente, para cada tipo de rede, o factor 2 
(alterável no decorrer do cálculo), obteve-se uma 
convergência bastante rápida, com uma redução 
na proporção aproximada de 1:0,7 do maior dos 
módulos das diferenças entre valores AU, em 
duas iterações consecutivas. 


Nalguns casos, no entanto, não se conseguiu 
convergência rápida, ou mesmo convergência, com 
um único valor de 2, sendo necessários pro- 
cessos mais elaborados e menos cómodos, por 
exemplo a consideração de valores diferentes 
para os vários nós, ou a utilização de processos 
de eliminação, ou a resolução directa do sistema 
de equações. 

Nalguns casos particulares, para alguns ele- 
mentos do circuito é possível definir um tensor 
impedância incremental, mas não é possível de- 
finir um tensor admitância incremental (por o 
determinante do tensor impedância incremental 
ser nulo). Nestes casos é possível a resolução do 
circuito a partir de equações de admitância, con- 
siderando um nó fictício adicional (ou deslocando 
eléctricamente, de maneira fictícia, um dos nós 
do circuito), de maneira a considerar admitâncias 
de ligação entre todos os nós figurados no cál- 
culo. 

Nos métodos de eliminação, utilizou-se a ge- 
neralização dos processos normais de eliminação 
de nós em circuitos definidos por impedâncias 
ou admitâncias complexas. 

Nos métodos de resolução directa de sistemas 
de equações, utilizou-se um processo de elimi- 
nação sucessiva com escolha da variável a eli- 
minar. 


4. EXEMPLO DE APLICAÇÃO AO CÁLCULO 
DE VARIAÇÕES DE TENSÃO NUMA REDE 


Como ilustração da aplicação dos tensores de 
comportamento incremental, considerou se um 
exemplo relativo ao cálculo de perturbações de 
tensão numa rede. Este exemplo é semelhante z 
um caso concreto estudado, mas foi propositada- 
mente simplificado para fazer ressaltar alguns as- 
pectos essenciais do comportamento da rede. A 
topologia base deste exemplo é definida na fig. 3. 

Essencialmente, pretende-se analisar o efeito 
nas tensões de correntes de perturbação numa 
zona industrial, com cargas sujeitas a variações 
bruscas importantes. 

Consideraram-se correntes de perturbações su- 
postas injectadas nos nós 4 e 5 e analisaram-se 
as perturbações de tensão correspondentes, na 
direcção da tensão anterior e na direcção orto- 
gonal. Fizeram-se várias hipóteses sobre o com- 
portamento incremental dos elementos da rede, 
no sentido de averiguar a influência dessas hipó- 
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teses no resultado, e em que medida e com que 
rigor é necessário conhecer e representar o com- 
portamento dos elementos da rede. 


() — jogo de barras 
— rede geral, suposta representável por 
impedâncias lineares e f.e.m. cons- 
tantes 
Baf — rede particular, tendo sido feitas sobre 
o seu comportamento as hipóteses in- 
dicadas 
— gerador 


— elemento de ligação entre nós da rede, 
do tipo de linhas ou transformadores 


—- elemento em derivação num nó da 
rede, do tipo de bateria de condensa- 
dores ou capacidade equivalente de 
linhas. 


Fig. 3 


A título ilustrativo, apresenta-se, nos quadros 
1 a 4 um extracto dos resultados em que se in- 
dicam, em unidades base (!), a variação de ten- 
são em diversos pontos, no sentido da tensão 
anterior nesses pontos (AU,) e em sentido orto- 
gonal (AUq), para perturbações de corrente de 
derterminada amplitude no ponto 4 (quadros 1 
e 2) ou no ponto 5 (quadros 3 e 4),na direcção 


= 


(!) Interessa sómente, para o efeito deste artigo, a 
comparação entre pontos ou hipóteses diferentes e não o 
valor absoluto das grandezas. 
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da tensão no ponto onde se dá a perturbação de 
corrente (quadros 1 e 3) e em direcção ortogonal 
(quadros 2 e 4). 


Consideram-se em alternativa: 
a) para os circuitos B, C, D, E e E: 


a.1) representação por admitância incremental 
igual a F vezes a admitância média, de- 
finida pelo cociente (complexo) da cor- 
rente pela tensão (em representação com- 
plexa) em condições base, tendo F, em 
alternativa, os valores 1; 0; — 0,5; — 1; 


a.2) representação por tensor admitância ou 
impedância incremental correspondente a 
elasticidade de potência activa £: e de po- 
tência reactiva n em relação ao módulo 
da tensão, supondo : == ne em alternativa 
igual a 2; 1; 0,5; 0; 


b) para o circuito G: 


b.1) comportamento análogo ao dos circuitos 
Ba E: 


b.2) comportamento correspondente a potência 
activa e módulo da tensão constante (hi- 
pótese representativa de um gerador com 
regulação de tensão e de potência a valores 
constantes). 


Tem-se assim, para cada quadro, um conjunto 
de 16 hipóteses, referenciadas de 1 a 16. 

A comparação dos valores situados numa 
mesma coluna de cada quadro indica a diferença 
entre os resultados para as diferentes hipóteses. 

Convém notar que a hipótese 12 corresponde, 
mercê das condições particulares dos parâmetros 
e da topologia da rede, à vizinhança de uma 
situação tipo «ressonância», o que explica os gran- 
des valores de perturbação, o mesmo sucedendo, 
embora menos acentuadamente, na hipótese 16. 

Verifica-se que, mesmo abstraindo das hipó- 
teses 12 e 16, as caracteristicas particulares dos 
elementos da rede e o método de representação 
têm enorme influência no resultado. 

Sem pretender estabelecer regras universais 
num domínio em que as situações podem variar 
muito de um caso a outro, pode concluir-se que: 


a) Naqueles casos, que podem considerar-se 
normais, em que as cargas tenham valores 
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b) 


c) 
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de : e 1 apreciávelmente diferentes de 2, e 
ou em que nas cargas ou complexos indus- 
triais haja proporção apreciável de má- 
quinas rotativas, é em geral necessário, 
para obter resultados válidos no cálculo de 
quedas de tensão correspondentes a per- 
turbação de carga, conhecer e representar 
com alguma aproximação o comportamento 
incremental das cargas e máquinas e a 
actuação dos elementos de regulação. 


A representação de elementos do circuito 
por admitâncias ou impendâncias incremen- 
tais (em sentido complexo) mesmo que 
corrigidas por um factor F, não é em geral 
e sob vários aspectos, satisfatóriamente 
correcta. 


Na medida em que, para efeitos de com- 
portamento incremental e variações de ten- 
são no sentido da tensão anterior, os dife- 
rentes elementos dominantes tenham carác- 
ter aproximadamente indutivo, isto é, os 
acréscimos de tensão e corrente nesses ele- 


mentos estejam aproximadamente em qua- 


dratura, as variações de tensão na direcção 
da tensão anterior, associadas a variações 
de corrente em quadratura com essa tensão, 
podem ser estimadas, em primeira aproxi- 
mação, representando os elementos com 
elasticidade : da potência activa e 7 da po- 
tência reactiva em relação ao módulo da 
tensão por admitâncias incrementais (em 
sentido complexo) de componentes activa 
e reactiva iguais aos valores correspon- 
dentes à situação base multiplicados respec- 
tivamente por:— 1 en—l, 


d) 


É difícil, no entanto, estabelecer «a priori», 
no caso de redes reais relativamente com- 
plexas, em que medida esta hipótese pode 
ser aplicada e prever a ordem de grandeza 
do erro cometido. 

A utilização dos tensores admitância ou 
impedância incremental permite, com gran- 
de facilidade, representar o comportamento 
correcto das redes, para pequenas pertur- 
bações, 


5. CONCLUSÃO 


Dos trabalhos efectuados no domínio dos 
assuntos tratados neste artigo, conclui-se desig- 
nadamente que: 


a) 


b) 


Os conceitos de tensor impedância e admi- 
tância incrementais são bastante cómodos 
e úteis no estudo de problemas de redes 
eléctricas, por exemplo no estudo de ques- 
tões de variações de tensão, variações de 
cargas, de perdas e de comportamento de 
redes tendo em conta o funcionamento 
de elementos de regulação. Estes conceitos 
constituem uma forma bastante simples de 
ter em conta efeitos e actuações que, dou- 
tra forma, na representação global de uma 
rede, é difícil ou laborioso considerar de 
maneira satisfatória. 

A representação tensorial de elementos de 
rede e os programas correspondentes ela- 
borados, quer para aplicação e estudos de 
comportamento incremental quer para ou- 
tras aplicações, afiguram-se bastante úteis 
e cômodos no estudo de redes eléctricas. 
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CONSTRUÇÃO CIVIL 


RESUMO 


Considerações sobre os processos tradicionais de 
construção civil e sua evolução para se adaptarem às 
diferentes condições econômicas. Novos processos, den- 
tro da técnica industrializada, e revisão dos metodos 
de trabalho, de organização das empresas de constru- 
(ão, e da própria concepção dos projectos, para se poder 
corresponder às necessidades da industrialização, 


1— Por «construção civil» designamos a acti- 
vidade humana relacionada com a construção de 
edifícios, sejam destinados a habitação ou a acti- 
vidades colectivas como hotéis, escolas, hospi- 
tais, etc. 

A construção civil tem mantido, na sua gene- 
ralidade, uma actuação de tipo «artezanal», uma 
vez que trabalha por encomenda, executa sempre 
modelos diferentes utilizando muita mão-de-obra 
e recorrendo pouco à energia mecânica. Várias 
razões concorreram, no passado, para essa si- 
tuação, o carácter regional das organizações li- 
gadas ao emprego dos materiais existentes no 
local, com utilização de mão-de-obra abundante, 
que fazia a sua aprendizagem segundo métodos 
tradicionais que exigiam uma longa prática e que 
eram transmitidos de geração em geração; os 
projectistas, trabalhando isoladamente, sem con- 
tinuidade e sem possibilidade de beneficiar da 
divisão do trabalho; os industriais desinteres- 
sados em aperfeiçoar as técnicas que utilizavam 
e, maturalmente, receosos em tentar aplicações 
vultuosas de capitais, face às incertezas dos mer- 
cados ; os serviços públicos, actuando muitas ve- 
zes com mentalidade ultrapassada e com falta de 
dinamismo; a opinião pública, particularmente 
respeitosa pelos processos já consagrados, tole- 
rante em relação aos defeitos, e defensora das 
virtudes da construção civil tradicional, etc. 

Nesta fase a actuação dos projectistas era 
muito reduzida, limitava-se aos elementos gerais 
do projecto e a indicações que iam dando aos 
construtores, à medida das necessidades. 
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INDUSTRIALIZADA 
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SYNOPSIS 


Considerations on the traditional building procedure 
and its evolution as to adapt to the different economic 
conditions. 

New procedure, inside the industrialized technique, 
and revision of working methods, building concerns 
organisation, and the conception of projects, as to cor- 
respond to the industrialization requirements. 


2— A situação atrás referida foi modificada 
com a industrialização dos transportes (ferro- 
viários e por estrada), que tornou técnica e eco- 
nômicamente possível o intercâmbio de materiais 
de construção civil entre as diversas regiões, des- 
truindo o equilíbrio de interesses anteriormente 
existente. 

Isso tornou possível, nomeadamente, a gene- 
ralização do emprego do cimento e a utilização, 
cada vez maior, do betão armado. Esta utilização 
foi também o resultado da criação de métodos 
industrializados para o fabrico, transporte e colo- 
cação deste último material, que permitiram 
passar a alcançar, sistemáticamente, e empre- 
gando pouca mão-de-obra, prazos de execução 
e preços inferiores aos anteriormente consegui- 
dos com a utilização da alvenaria de pedra. Com 
efeito, a extracção, o transporte e a colocação 
deste material em obra, exige larga aplicação de 
mão-de-obra especializada. 

Tornou ainda possível o desenvolvimento 
da préfabricação ligeira, que pôs à disposição 
dos construtores certos elementos de série como 
a estrutura de pavimentos de edifícios, asnas e 
cascas para coberturas, serralharias diversas, por- 
tas e janelas de determinados tipos, etc. 

Esta situação, aliada à maior complexidade 
das obras devido ao maior número de instala- 
ções de diversa natureza e ao maior número de 
materiais aplicados transformou a organização 
tradicional da construção civil. Esta passou, de 
uma instituição fortemente centralizada, que tudo 
executava, para uma organização descentralizada, 
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constituída por um número cada vez maior de 
gabinetes de estudo especializados e de sub-em- 
preiteiros independentes. O «arquitecto» e o «en- 
genheiro» encarregam-se da coordenação dos 
diversos trabalhos de projecto e o «engenheiro» 
ou o «construtor» da coordenação dos diversos 
intervenientes na construção dos edifícios, isto 
é, fornecedores de materiais e montadores, sub- 
-empreiteiros de mão-de-obra, e de trabalhos 
complementares. 

A resolução dos problemas postos nesta base 
exige, maturalmente, a existência de serviços 
complexos. A eficiência dos gabinetes de estudos 
e das empresas de construção é medida pela 
maneira como conseguem coordenar, de forma 
permanente, as diversas organizações que con- 
tribuem para a realização das obras. Isso exige 
já o recurso a um planeamento mais cuidadoso, 
em moldes novos. 


3 — À época actual põe problemas que antes 
não existiam, resultantes do grande aumento da 
população do Mundo (duplicação em cerca de 
40 anos) e da obrigação social de resolver os 
problemas materiais básicos de todos os homens, 
pondo à sua disposição os bens indispensáveis 
a uma via digna. 

A maneira de conseguir tal objectivo é atra- 
vés da industrialização, que tem vindo a abran- 
ger cada vez mais actividades humanas, sendo a 
construção civil o último baluarte a conquistar. 

Por «construção civil industrializada» enten- 
de-se a construção civil que utiliza métodos de 
técnica industrial com vista a estender, a esta 
actividade, o objectivo da industrialização que 
se traduz, neste caso, em obter edifícios em 
grande quantidade, a menor prazo, a preço mais 
baixo e de melhor qualidade. 

Para conseguir tais resultados é indispensável, 
antes de mais nada, ter assegurado um mercado 
potencial suficiente para permitir à organização 
trabalhar com continuidade. Com efeito, não é 
possível amortizar em poucas unidades o custo 
das instalações industriais e não é viável montar 
em boas condições uma organização especiali- 
zada sem garantias de continuidade. Consi- 
dera-se que só é viável montar uma indústria de 
construção civil com a certeza de um mercado 
mínimo inicial de 4000 a 5000 fogos executando 
anualmente cerca de 1000 fogos, e repetindo-se 
cada tipo pelo menos 100 a 150 vezes. 
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Esta questão prende-se, naturalmente, com à 
política nacional de investimentos porque não é, 
em regra, possível, deixar à iniciativa privada o 
financiamento global necessário à resolução do 
problema da habitação, dados os montantes a 
dispender e a prática impossibilidade de reem- 
bolso a curto prazo, contando só com os directos 
interessados. 

A produção industrial efectua-se numa fábrica, 
eliminando o trabalho sob a acção das intempé- 
ries e reduzindo a discontinuidade de emprego 
e a deslocação ou nomadismo, de um estaleiro 
para o outro, que caracteriza a construção tradi- 
cional. 

As dimensões de um edifício são assaz gran- 
des, da ordem de uma a muitas dezenas de me- 
tros, portanto muito superiores ao que as fábri- 
cas costumam produzir. Há, deste modo, que 
sub-dividir o edifício em elementos mais peque- 
nos, cujo transporte seja viável, e cuja ligação 
se possa processar, na obra, em boas condições. 
Com efeito, dado o peso e o volume de uma 
construção, não é possível pensar na saída dos 
edifícios da fábrica já depois de completamente 
montados. 

A produção industrial implica uma certa «repe- 
tição», isto é, não pode considerar a realização 
de um modelo único, embora possa executar 
qualquer modelo desde que ele se repita um nú- 
mero interessante de vezes. Convém aqui salien- 
tar que a repetição que se refere diz respeito aos 
«elementos» de construção e não, necessaria- 
mente, ao conjunto dos edifícios. 

Com efeito, a introdução do trabalho das má- 
quinas só é rendável na produção em série, isto é, 
na produção de produtos iguais durante um pe- 
ríodo de razoável duração. É geralmente consi- 
derado que, nos trabalhos artezanais da cons- 
trução, se obtém uma redução de custo de 10 a 
20 “/o quando se repete dez vezes um mesmo 
trabalho. Nas actividades industriais este aspecto 
é, evidentemente, mais relevante, pois que não 
é possível uma organização económica do tra- 
balho das máquinas e do pessoal se as tarefas a 
executar variarem constantemente. 

Mas não é só na organização do trabalho que 
a repetição tem vantagens. É possível dispender 
mais dinheiro no projecto de uma grande uni- 
dade — levando muito mais longe o estudo — do 
que no caso de pequenas unidades que se não 
repetem. Dai um muito maior incentivo para a 
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